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Control of electron spin decoherence in contact with a mesoscopic bath of many interacting nuclear 
spins in an In As quantum dot is studied by solving the coupled quantum dynamics. The nuclear 
spin bath, because of its bifurcated evolution predicated on the electron spin up or down state, 
measures the which-state information of the electron spin and hence diminishes its coherence. The 
many-body dynamics of nuclear spin bath is solved with a pair-correlation approximation. In the 
relevant timescale, nuclear pair-wise flip-flops, as elementary excitations in the mesoscopic bath, can 
be mapped into the precession of non-interacting pseudo-spins. Such mapping provides a geometrical 
picture for understanding the decoherence and for devising control schemes. A close examination of 
nuclear bath dynamics reveals a wealth of phenomena and new possibilities of controlling the electron 
spin decoherence. For example, when the electron spin is flipped by a 7r-pulse at r, its coherence will 
partially recover at \/2t as a consequence of quantum disentanglement from the mesoscopic bath. In 
contrast to the re-focusing of inhomogeneously broadened phases by conventional spin-echoes, the 
disentanglement is realized through shepherding quantum evolution of the bath state via control of 
the quantum object. A concatenated construction of pulse sequences can eliminate the decoherence 
with arbitrary accuracy, with the nuclear-nuclear spin interaction strength acting as the controlling 
small parameter. 

PACS numbers: 03.65.Yz,03.67.Pp,76.30.-v, 71.70.Jp 



I. INTRODUCTION 

A quantum system, unlike a classical one, can be 
in a coherent superposition of constituent states. This 
coherence is the wellspring of quantum properties and 
key to quantum technology. The contact of a quan- 
tum object witha macroscopic system causes loss of 
state coherence.'^^^El Advances towards quantum technol- 
ogy have effectively substituted the macroscopic environ- 
ments by mesoscopic ones. Usually, the contact inter- 
action between the quantum object and a "particle" of 
the bath weakens with increasing the bath size (defined 
by the number of particles, N), while the interaction be- 
tween particles inside the bath is independent of the bath 
size. When the bath size is in the mesoscopic regime such 
that the object-bath interaction dominates their interac- 
tion with the rest of universe, the quantum object and 
the mesoscopic bath evolves as a closed system in the rel- 
evant timescales. Thus, a proper description of the quan- 
tum object in such a mesoscopic bath is given by the full 
quantum mechanical solution of the coupled object-bath 
evolution,!^ as opposed to the semiclassical treatment in 
conventional decoherence studies.!^ 

Spins of single electrons confined in semiconductor 
quant um dots are paradigmatic systems in me soscopic 
physicJSS and in spin-based quantum technologjff^IIIIill^ 
Spin decoherence is a main limiting factor to quan- 
tum properties and has been extensively studied both 

in thcoric d'^'^'^^'^'^'^^'^^'^'^'^^'^^'^"'^^'^^^'^^'^'^'^^'^^'^'^'^ 
and in experiments.'22l 3" l 3i | 32 | 33 | 34 | 35 | 36 | 37 | 38 | 39 | 

has been well estab Ushed by theorieJ^^E^OSIlI] and 
experiment d^^|30|3i | 32 | ^]-^g electron spin decoherence 



caused by phonon scattering is negligible at a tempera- 
ture lower than a few Kelvins. The relevant bath at low 
temperature then is a spin bath^*^ composed of lattice 
nuclear spins in a quantum dot (QD). The electron spin 
is coupled to a nuclear spin through the contact hyper- 
fine interaction with magnitude inversely proportional 
to the total number N of nuclei in the QD. For the QD 
size of interest, the hyperfine coupling is much stronger 
than the mutual interaction between nuclear spins. 
Therefore, a mesoscopic bath consisting of all nuclear 
spins within the QD (i.e., in direct contact with the 
electron) is identified. Electron spin decoherence at low 
temperatures is determined by the quantum dynamics 
of the mesoscopic spin system. 

Under a moderate to strong magnetic field {B > 
0.1 T), the Zeeman energy of the electron is by orders 
of magnitude larger than the hyperfine coupling, so that 
the longitudinal electron spin relaxation by off-diagonal 
hyperfine interaction is virtually suppressed.** The Hamil- 
tonian of the electron-nuclear spin system can then be 
reduced to the form H = \±)H±{±\ where |±) de- 
note the electron spin eigenstates in the external field and 
H± are the nuclear bath Hamiltonians depending on the 
electron spin states. The pure dephasing (decoherence) 
of the electron spin is caused by the electron-nuclear spin 
entanglemenlj^ as follows. Let the electron-nuclear spin 
system start from a product (i.e. unentangled) state, 
(C+I+) + C_|-)) «) \ J). The nuclear spin state \J) 
would be driven by the Hamiltonians H± to the states 
iJ'^it)) corresponding to the electron states |±), respec- 
tively. Thus the electron-nuclear spins would evolve into 
an entangled state, C+\+) «) \ J+{t)) + C-\-) (g) \J~{t)). 
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FIG. 1: Schematic illustration of the bifurcation of the bath 
state evolution conditioned on the electron spin up or down 
state and the exchange of the evolution direction of the bath 
pathways when the electron spin is flipped. 

The ofF-diagonal element of the reduced density matrix 
of the electron spin = C* C+ {J-{t)\J+{t)) mea- 

sures the electron spin coherence. Therefore, the bifur- 
cated evolution of the nuclear bath leads to the loss of 
electron spin coherence. From the viewpoint of quantum 
measurement, the electron spin states are registered by 
different nuclear bath states and a measured quantu m ob- 
ject sustains no coherence in the measurement basis .l^^IiH 
Within the timescale when the irreversible leakage of 
the quantum coherence of a mesoscopic bath into the 
macroscopic environment is negligible, it is still possi- 
ble to control the mesoscopic bath dynamics'^ and hence 
the decoherence of the quantum object embedded in the 
bath.l^ As illustrated in Fig. [l] the bath evolution path- 
ways in the Hilbert space corresponding to opposite elec- 
tron spin states will exchange their evolution directions 
when the electron spin is flipped: 

{C+\+) + C^\-))®\J) 
C+\+)(g>\J+it)) + C-\-)(g>\J-{t)) 
C+\-)(^\J+it')) + C-\+)(^\J~{t')). 

Thereafter, the two bath pathways \ J'^{t')) and \ J^{f)) 
intersect at some later time, i.e., \J^~^{t')) = \J^~{t')) — 
\J^{t')). At the intersection, the electron spin is disen- 
tangled from the bath and its lost coherence is recovered 
(by the controlled erasure of the quantum information 
registered in the bath). Such a discovery of the rever- 
sal of coherent dynamics can be traced back to the early 
theoretical a nd exp erimental findings on the Loschmidt 
echo in NMR.EsEa 

The recovery of the lost coherence (recoherence) 
by quantum disentanglement is fundamentally different 
from the conventional spin echo realized by the refocusing 
of the random phase in an inhomogeneously broadened 
ensemble. In general, the disentanglement can occur at 
a time different from the spin echo time. For example, 
by a single flip of the electron spin at t = r, there will 
be a prominent recoherence at a magic time as a 
consequence of disentanglement^ while the conventional 
spin echo would occur at 2r. When the electron spin 



is observed in ensemble measurement (realized either by 
using many similar QDs or by cycling measurements on 
a single dot), the disentanglement-induced recoherence 
will be concealed by the inhomogeneous broadening un- 
less it is forced to take place at a spin echo time by proper 
design of pulse sequences. 

The recoherence by disentanglement can be ex- 
ploited as a coherence protection scheme in lieu 
of the dynamical decoupl ing schemes developed for 
NMR spectr oscopiej^* *^ ! ^^ ! and recently for quantum 
computation.EHIinilllllllllll While the dynamical decou- 
pling schemes seek to eliminate the object-bath inter- 
action through rapid rotation of a quantum object, the 
disentanglement focuses on controlling the wavefunction 
evolution of the bath and in general does not lead to a 
vanishing object-bath coupling. It will be shown that 
in the disentanglement scheme, the controlling small pa- 
rameter for coherence protection is determined by the 
interactions within the bath instead of the object-bath 
coupling, in contrast with the dynamical decoupling 
schemes. 

This paper is organized as follow: Sec. |TT] describes a 
speciflc model system (a self-assembled InAs QD) and 
summarizes the theoretical ingredients for solving the 
electron-nuclear spin dynamics. The free-induction decay 
(FID) in single-system dynamics (without ensemble aver- 
age) and the Hahn echo in ensemble dynamics are studied 
in Sec. |III| and Sec. |IV| respectively. Sec. [V] presents a 
close examination of the single-system dynamics in the 
single-pulse Hahn-echo configuration, revealing a coher- 
ence recovery due to disentanglement at a magic time 
which would otherwise be invisible in ensemble-averaged 
signals. Sec. |VI| studies the disentanglement under con- 
trol of pulse sequences, in particular, the concatenated 
sequences. Further discussions including the comparison 
between the quantum disentanglement and the dynam- 
ical decoupling are presented in the summary Section. 
Some technical details are given in the Appendices. 

II. THEORY 

A. The Model 

The system consists of an electron with spin-1 /2 Sc and 
N nuclear spins, 3n,a, with Zeeman energies f2o and Ua 
under a magnetic field iJoxt , respectively, where n denotes 
the position and a denotes the isotope type. Hereafter 
the subscript a will be absorbed into n unless it is needed 
for clarity. InAs has the Zincblende structure, with In 
and As ions located in two interpenetrating face-centered 
cubic lattices. The natural isotope abundance in InAs 
materials is 100%, 4.28% and 95.72% for ^^As, ^^^In, 
and ^^^In, respectively. All the isotopes have non-zero 
nuclear spin moments, namely, Jtsas = 3/2 and JnHn — 
Jii5in = 9/2. The self-assembled InAs QD under typical 
growth condition is modeled as a rectangular quantum 
box with the growth direction along [001] and the in- 
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TABLE I: InAs material parameters, where 7^ is the nuclear 
gyromagnetic ratios, 7™ the effective nuclear gyromagnetic 
ratios for the pseudo-exchange interaction, 7e the electron gy- 
romagnetic ratio, 7I the effective electron gyromagnetic ratio 
in the QD, and Ac the hyperfine constants. 





'^As 


^^^In 


^^^In 


Abundance 


100% 


4.28% 


95.72% 


Spin moment Ja 


3/2 


9/2 


9/2 


7a (10« s-^T-^) 


-F45.8 


-)-58.5 


4-58.6 


7^" (10'' s-^T-^)" 


-H34.0 


-1-70.1 


4-70.2 


Ac. (10^ s-^)" 


-f69.8 


-f85.1 


4-85.3 


7, = +0.176, 7e* = -0.132 (10"^ s~"T'")" 


lattice constant cq = 6.06 A 



"See details in Ref.|l[56] andlSTl 

'Estimated with the method in Ref. 1581 
'^g-factor taken to be —1.5. 



TABLE 11: The characteristic energy scales in an InAs QD 
with dimensions 35 x 35 x 6 nm (N ~ 0.3 X 10**), under a 
field of Bcxt ~ 10 Tesla. 



Units" 






En 








10" 


500 


1 


10-* 


10"" 


AteV 


10^ 


0.5 


10"^ 


10-'^ 


lo-'^ 


mK 


10* 


5 


10"^ 


10"'' 


lO"** 



"Throughout this paper, the units are chosen such that the Plank 
constant h and the Boltzmann constant ks are unity. 



FIG. 2: (a) Schematics of a quantum box model for an InAs 
QD. (b) Schematics of an electron (the shadow) and one layer 
of nuclear spins in the QD. The two boxes in dotted lines indi- 
cate two possible choices of boundary of the nuclear spin bath, 
which are relatively arbitrary due to the interaction between 
nuclei within and without the boundary. When the hyper- 
fine interaction dominates over the nuclear spin interaction, 
such arbitrariness has negligible effects on calculation of the 
electron spin decoherence as long as all the nuclei in direct 
contact with the electron spin have been enclosed. 

plane extension directions [110] and [lIOj.'^The electron 
is assumed to be confined by hardwall potential and the 
envelope wavefunction of the ground state is 



2 Trr^ 

La La 



0(ia-2|ra|): 



(1) 



Virtual flip-flops of the electron spin, however, could me- 
diate an extrinsic interaction between nuclear spins even 
when they are well separated in space (as illustrated in 
Fig. [3|.^ Furthermore, the non-secular part of the nuclear 
spin interaction which does not conserve the total Zee- 
man energy, including flip-flops of hetero-nuclear pairs 
and single-spin flips, can be neglected. Thus the total 
effective Hamiltonian is reduced to the form 



H — He + _ff ^ 



(±1 



(2) 



in the limit of long longitudinal relaxation time (Ti — 
cx)), where |±) are the eigenstates of S^, He = ^cS^ 
= ^nJnj the nuclear spin interaction 



where ra is the coordinate in the a direction([110], [110], 
or [001]), and L^ is the dimension of the dot along the 
indicated direction. The model system is illustrated in 
Fig.|2] 

The first-principles Hamiltonian for the electron- 
nuclear spin system includes the electron-nuclear hyper- 
fine interaction and various intrinsic nuclear-nuclear in- 
teractions. A detailed discussion was given in Ref. S) 
Under a moderate magnetic field (Bcxt > 0.1 T), the elec- 
tron spin fiip is virtually suppressed due to the large Zee- 
man energy mismatch between the electron and nuclei. 



H± = ±Ha + Hb + Hd±He 



(3) 



FIG. 3: Hyperfine-mediated interaction between two distant 
nuclear spins via the electron (top line). 
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with 



Ha 
Hb 
Hd 
He 



V ' "-"""^ f+ J- - V ,4 f+ J- (4-a^ 



n<m 



(4b) 
(4c) 
(4d) 



Here the summation with a prime runs over only the 
homo- nuclear pairs, Ha denotes the extrinsic interac- 
tion mediated by the hyperfine interaction, Hb denotes 
the ofF-diagonal part of the intrinsic nuclear interac- 
tion and Hd the diagonal part, and He is the diagonal 
part of the contact hyperfine interaction with amplitude 
On = -4aCo |/(R-n)|^- The material parameters are given 
in Table IT] and the typical energy scales are estimated 
in Table In] The hyperfine interaction has a typical en- 
ergy scale En ^ a„ ~ 10^ s^^ for a dot with about 10^ 
nuclei. The intrinsic nuclear spin-spin interaction is ef- 
fectively finite-ranged with the near-neighbor coupling 
Bn.7n ^ Dn.m ^ 10^ s^^. The cxtriusic interaction is 
"infinite-ranged" (coupling any two nuclear spins within 
the QD) and has opposite signs for opposite electron spin 
states. The extrinsic interaction, in contrast with the in- 
trinsic one, depends on the external field, and has an 
energy scale An,m ~ 1 to 10 s~^ for field strength i?ext 
varying from 10 to 1 Teslas. 



Due to the block diagonal form of the reduced Hamilto- 
nian Eq. (J2| , the correlation function has following prop- 
erties 



e iiN 



(9a) 
(9b) 

(9c) 



where the last equation has been expressed for the free- 
induction decay and can be straightforwardly extended 
to the dynamics under pulse control. No longitudinal re- 
laxation remains (Ti = oo) after the elimination of the 
electron spin fiip processes by the magnetic field. Since 
the Zeeman energy Hq results in only a global phase- 
factor e~*^=*, from now on we will consider only the in- 
teraction Hamiltonian Hi = H — Hq by dropping the 
trivial phase-factor. 

Since the nuclear Zeeman energy under the external 
field dominates over the hyperfine interaction and the 
nuclear spin interaction and the temperature to be con- 
sidered is much higher than the interaction energy, the 
thermal nuclear spin state can be taken as 



J 



(10) 



where IJ) = |j„) is an eigenstate of H^ and Pj 
YinPjn is the probability distribution with 



B. General formalism 

The electron-nuclear spin system is assumed initially 
prepared in a product state,^ 



(5) 



where the nuclear spins are in a thermal state with tem- 
perature T. The electron spin can start from a pure state 
or a mixed state.l^The state evolution is determined by 
the Louville equation 



dtm = AH,m]. 



(6) 



The reduced density matrix of the electron spin is ob- 
tained by partial trace over the nuclear spins as 



p%t) = TrN/5(t), 



(7) 



which is related to the initial state through the correla- 
tion superoperator C, namely 



L,(t)- ^/:^,.;^.,.,(i)p°M',.'(o), 



(8) 



where p^^ = (mIpIj^), and |/^), e {|+), |-)}. The 
electron spin relaxation is quantified by the correlation 
function C^^y-^iy. 



.IT 



(11) 



giving the population of the single-spin state |j„). We 
note that such a nuclear spin state has no off-diagonal co- 
herence (between the Zeeman energy eigenstates) . The 
correlation function of the electron spin coherence, ex- 
pressed in terms of the wavefunction overlap of the nu- 
clear spins, is 



{t)\J^{t)) 



(12) 



with \ J^{t)) = exp {^—iH±tj \ J)- The correlation func- 
tion £-1-, (^) is c losely related to the Loschmidt echo in 
literature.lSiE^I^ Being independent of the electron spin 
initial state, £+._(i) will be equated with the electron 
spin coherence without causing confusion. The deco- 
herence by entanglement is transparent from Eq. (12 1: 
With the nuclear spin states \J^{t)) driven by differ- 
ent Hamiltonians H±, the electron-nuclear spin state 
C+I+) (g) \ J^{t)) + C-\-) (g) \ J~{t)) becomes an entan- 
gled one when {J~{t)\ J^{t)) < 1, so the electron spin 
coherence is lost. 
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C. Essential assumptions 



D. Ensemble and single-system dynamics 



The method to be used in this paper consists of two 
steps. First, the electron spin decoherence in ensemble 
dynamics is factorized into two parts: the dephasing due 
to static inhoniogcneous broadening, and the decoherence 
due to the dynamical entanglement between the electron 
and the nuclei in single-system dynamics. Second, the 
nuclear spin dynamics is solved with the pair-correlation 
approximation. The following assumptions are essential 
for the theory: 

1. Mesoscopia: First, the QD should be small enough 
to guarantee the dominance of the hyperfine inter- 
action over the nuclear spin interaction so that the 
isolation of a mesoscopic bath is possible. Second, 
the bath size (i.e., the number of nuclear spins in 
the QD, N) should be large enough to have genuine 
decoherence (i.e., the Poincare time is much greater 
than the spin relaxation time Ti and T2). Third, 
N ^ y/N is required so that the central limit theo- 
rem in statistics can be utilized for the factorization 
of the single-system dynamics from the inhomoge- 
neous broadening effect. And finally, to justify the 
pair-correlation approximation, the bath size N is 
bounded so that the number of nuclear pair-flips 
in the relevant timescale is small compared with N 
(see Appendix [B|. 

2. Moderate to strong field: The electron Zeeman en- 
ergy f2e be much greater than the hyperfine e nergy 
a„ so that the electron spin flip is suppressed. EES 

3. Finite temperature: The temperature be higher 
than or comparable to the nuclear Zeeman energy 



(T 



> 



1 mK). Otherwise, in a near fully po- 



larized nuclear spin bath under an extremely low- 
temperature, the pair flip-flops occur only among 
a few available nuclear spins and hence are highly 
correlated, in which the pair-correlation approxi- 
mation is invalid. Also, the temperature should be 
low enough to prevent the phonon scattering effect 
(r < 1 K) 

4. Short-time dynamics: The timescales under consid- 
eration should be short compared with the inverse 
nuclear spin interaction strength (~ 1 ms) so that 
the number of pair-flips over time is small compared 
with N. The timescales of interest in the decoher- 
ence problem include the coherence memory time, 
pulse delay time in control sequences, and the total 
duration of a control pulse sequence. 

The other approximations in the model, including the as- 
sumptions on the QD shape, the electron wavefunction, 
and the specific form of the spin interactions, are incon- 
sequential to the theory. 



The electron spin decoherence given in Eq. ^ includes 
two contributions, namely, the thermal fluctuations due 
to the ensemble average of different nuclear spin configu- 
rations and the quantum entanglement due to 
the dynamical evolution starting from a pure state \J'), 
which is conditioned on the electron spin state. We use 
the term single-system dynamics to denote the quantum 
evolution governed by the Schrodinger equation 



idt\J^it))^H±\J^{t)), 



(13) 



and ensemble dynamics to denote the ensemble average 
of the single-system dynamics. 

The ensemble dynamics was studied in the for- 
malism of density matrix with the pair-correlation 
approximation.'2Sl Many important features of the quan- 
tum entanglement process, however, could have been 
overlooked if without resolving the single-system dynam- 
ics since the static inhomogeneous broadening usually has 
much stronger effects on the decoherence than the dy- 
namical entanglement does. We choose a different pro- 
cedure to study the problem. Namely, we will first solve 
the single-system dynamics starting from a certain nu- 
clear spin configuration, and then construct the ensemble 
dynamics via 



Cl_{t)^{J-it)\J+it)), 



(14a) 
(14b) 



J 



Such a construction of the ensemble dynamics, in general, 
would require sampling a large number of initial states 
\J') in the thermal ensemble. 

When the system is sufficiently large and the tempera- 
ture is appreciable for the nuclear spins, the decoherence 
due to quantum fluctuations in single-system dynamics 
(Eq. 14a I is almost the same for all possible initial nu- 



clear configurations in the thermal ensemble, except for a 
global phase factor related to the static Overhauser field 
resulting from the diagonal hyperfine interaction. This 
is confirmed by numerical verification (see Appendix [A| . 
Here we give a physical argument: For a sufficiently large 
number of nuclear spins far from being fully polarized, 
the number of spins available for pair-wise flip-fiops is 
large, so the electron spin decoherence is virtually deter- 
mined by the nuclear spin excitation spectra (or density 
of states). By the statistical central limit theorem, the 
excitation spectrum of a large system is the same (up to 
a relative variance ~ l/\/iV) for different initial states 
IJ). Then the electron spin decoherence can be factor- 
ized into a contribution from the static inhomogeneous 
broadening and one from the quantum fiuctuation 
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in single-system dynamics defined as 



Cl^_{t) EE \{J-it)\J+{t))\ 



lE.jt 



J 



(15a) 
(15b) 

J P{E)e-'^'dE, (15c) 



where the hyperfine energy due to the local Overhauser 
field in a certain nuclear spin configuration is 



E 



J 



E 



Jn^n 1 



(16) 



and the inhomogeneous broadening distribution is de- 
fined as 



P{E) = Y,Pj5{E-Ej). 
J 



(17) 



The quantum fluctuation effect in the single-system dy- 
namics can be evaluated for an initial configuration \J') 
randomly selected from the ensemble, with the global 
phase factor absorbed into the inhomogeneous broaden- 
ing. 

The effect of the static inhomogeneous broadening can 
be calculated by neglecting the nuclear spin interaction 
(setting Ha — Hb — 0). Without the fiip-fiop dynam- 
ics, the nuclear spin ensemble would be a random dis- 
tribution of frozen configurations, which, as a standard 
multinomial distribution for a large system, leads to a 
Gaussian distribution of the Overhauser field 



PiEj) 



1 



2ttT* 



(18) 



where the averaged local field Eq and the inhomogeneous 



0.2 



0.1 



0.0 



£■„ = 0.345 I 




r;=i.76 










V 



-6 



-3 3 

£(10" s"^) 



FIG. 4: Distribution of the Overhauser field for a thermal en- 
semble of nuclear spins in an InAs QD of size 34 x 34 x 3 nm^ 
at a temperature of 1 K and under an external field of 1 Tesla. 
The dots in the figure are obtained by numerical simulation 
of 40000 random samples, and the curve is calculated with 
Eqs. ( |18|19[ ). The arrows indicate four randomly selected nu- 
clear spin configurations to be used for Fig.[l9]in Appendix]^ 



broadening are 



(19a) 



n = /5](47VaA2iV-2j2_ij2^Ar-i), (19b) 



respectively, with ja = J2jjPa,j, Ja = HjfPotj, and 
Na denoting the number of a- type nuclei {N = Y^^ Na). 
A typical example of the Overhauser field distribution is 
shown in Fig.|4] The electron spin dephasing by inhomo- 
geneous broadening is a Gaussian decay 



-iEot-{t/T;Y 



(20) 



with the effective dephasing time Tj* = "/^/'^2- The in- 
homogeneous broadening effect can be removed by spin 
echo. In the rotating frame rested on the electron spin, 
the Overhauser field changes its sign each time the elec- 
tron spin is fiipped by a short pulse. Thus under the fiip- 
ping operation of a sequence of pulses applied at t = ii, 
^2, ■ • ■, tn, the dephasing by the inhomogeneous broad- 
ening becomes 



J 



(21) 



with t = ti-{t2- h) + {h-t2) (-l)"(i - tn)- Ob- 
viously, the phase accumulation from the random Over- 
hauser field is cancelled at i = 0, leading to the spin echo. 
At the spin-echo time, the electron spin decoherence re- 
sults solely from the dynamical quantum entanglement. 

Under typical conditions considered in this paper, 
is found to be of the order of nanosecon ds, consist ent 
with various experimental measurements. '^^'^^'^lllll As 
will be shown bel ow, in agreement with the available 
experimental data,'^^'^SIMI the electron spin decoherence 
by quantum fluctuations in single-system dynamics has 
a timescale in the order of microseconds. Thus, from 



Eq. ( 15 1, the ensemble dynamics will be dominated by in- 
homogeneous broadening and the single-system dynam- 
ics is virtually invisible in ensemble experiments, unless 
the inhomogeneous broadening effect is removed by spin 
echo. The decay of spin-echo signals is usually attributed 
to the pure decoherence due to quantum fluctuations. We 
have shown, however, by examining directly the single- 
system dynamics, the flipping pulses have non-trivial ef- 
fects on the nuclear spin dynamics and, therefore, on 
the dynamical entanglement, making the spin-echo decay 
time significantly different from the FID time in single- 
system dynamics."* 

From now on, we will concentrate on the single-system 
dynamics, and show a few surprising effects which would 
otherwise be concealed by inhomogeneous broadening in 
ensemble dynamics. It should be emphasized that the 
single-system dynamics in a QD at a finite temperature 
is not a mathematical idealization but has measurable 
effects. For example, a projective measurement of the 
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local Overhauser field could be used to limit the nuclear 
spin configurations by post-selection and t hus to obs erve 
single-electron dynamics without spin echoJ^^^^^lSlIlSI The 
single-system dynamics is the basic to quantum technolo- 
gies such as quantum computation which can not be per- 
formed in ensembles for scalability to large systems with- 
out exponential explosion in resourceP^l 



tions. The transition for the pair-fiip driven by the op- 
erator J^J~ is 



\jn)n\jni)r 



\jn + 



(22) 



The transition between the electron-nuclear spin many- 
body states through the pair-fiip is denoted by 



E. Pair-correlation approximation and pseudo-spin 
model for nuclear spin dynamics 

With the Hamiltonian given in Eq. Q, the quantum 
evolution of the nuclear spins is driven by the pair-wise 
homo-nuclear spin fiip-fiops as the elementary excita- 



\±)\J)^\±)\J,k), 



(23) 



where k is shorthand for the pair-fiip in Eq. ( 22 1 . The 



transition described in Eq. (23 1 is characterized by the 



matrix elements ±Ak + Bk and the energy costs Dk^Ek, 
which are derived from the microscopic models as follows 



Ak = {J, k\ Ha \J) = -J^VJn {Jn + 1) - in (j„ + l)^/Jm {Jrn + 1) - jrn (jrn - 1), (24a) 

Bk = {J, k\ Hb \J) = B„,,„V^« {Jn + 1) - in Un + I) V Jm {J,n + 1) - jm (jm - 1), (24b) 
Dk = {JMHD\J,k)^{J\HD\J) 

^ ^ ^ n.7i' jn' ~^ ^ ^ -^m.rn'Jm' ^ ^n.m ~^ ^n^n {jn "h 1) ^ ^m,m (jm 1) ? (^^^) 

Ek = {J,k\HE\J,k)-(J\HE\J) = ia,,-a„,)/2. (24d) 



The basic processes for the electron spin decoherence 
may be described as follows: The off-diagonal part of 
the nuclear spin interaction (including the intrinsic one 
Hb and the extrinsic hyperfine-mediated one H^) causes 
the pair- wise nuclear spin fiip-fiops, leading to a fiuctuat- 
ing local Overhauser field and in turn a random dynam- 
ical phase for the electron spin. In the quantum picture, 
the entanglement is developed because the quantum evo- 
lution of the nuclear spins driven by the Hamiltonians 
±Ha + Hb + Hd ± He is conditioned on the electron 
spin state through the ± signs originating from the hy- 
perfine interaction. 

As established in Ref. i4| and discussed in further de- 
tails in Appendix [C] for pulse sequence controls, within a 
time t much smaller than the inverse nuclear interaction 
strength, the total number of pair-fiip excitations Nfnp is 
much smaller than the number of nuclei N. The prob- 
ability of having pair-fiips correlated is estimated to be 
^corr 1 — e~*^''"p/^ {q being the number of nearest 
neighbors), which, as also shown by a posteriori numeri- 
cal check (see Appendix [C| , is bounded by ~ 10% in the 
worst scenario studied in this paper. Thus, the pair-fiips 
as elementary excitations from the initial state can be 
treated as independent of each other, with a relative er- 
ror e < Pcorr (see Appendix [d| . Then the single-system 
dynamics \J'^{t)) can be described by the excitation of 
pair-correlations as non-interacting quasi-particles from 



the "vacuum" state \J'), driven by the "low-energy" ef- 
fective Hamiltonian, 

^l-E^^^Ex^^^A (25) 

k k 

which has been written in such a way that each pair-fiip 
is treated as a pseudo-spin 1/2, represented by the Pauli 
matrix o^j., with k labelling all possible pair-fiips. The 
time evolution from the initial state \J') can be viewed 
as the rotation of the pseudo-spins, initially all polarized 
along the —z pseudo-axis: 0^ | [)k, under the effective 
pseudo-magnetic field, 

Xk=i.±'^Ak + 2Bk,Q,Du±Eu), (26) 

for the electron spin state |±), respectively. Such a treat- 
ment of the nuclear spin correlations amounts to taking 
into account all the pairwise correlations and neglecting 
the higher-orders, as justified by a systematic study based 
on linked cluster expansion .^3 

The entanglement between the electron spin and the 
pseudo-spins, and hence the electron spin decoherence, 
are developed as the pseudo-spins precess about different 
pseudo-fields Xk corresponding to different electron spin 
states |±). Namely, the product state 

(C+|+) + C_h))(g)U),, (27) 

k 
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FIG. 5: The Bloch vectors for the pseudo-spin states 
precess about the pseudo-fields Xk f*^'' electron spin states 
±), respectively. 



will evolve into an entangled one 

C+\+) (g) \i4{t)) + CA-) (g) |^,-(t)) , (28) 

k k 

with 

IV'^ W) = exp {-^xtt) li) = e-^^i-^-'/'\ I). (29) 

The electron spin coherence is determined by the overlap 
between the "conjugate" states |'0fc^(i)) and |^jr(i)) of 
the pseudo-spins, 

^■+At) = l[\(^kit)\i^tit))\- (30) 



The pseudo-spin states have a geometrical representa- 
tion as Bloch vectors, 



.± — /„;.± 



(31) 



Fig. [5] shows the evolution of two conjugate pseudo-spin 
states IV'fe (^)) in two paths of the conjugate Bloch vec- 



tors cr^ on a sphere. The overlap between the conjugate 



pseudo-spin states is related to the distance between the 



two vectors Sh 



by 



|(^fe(0| V^,+ (i))| 



Si 



(32) 



Thus, the square of the distance 5| is the distinguishabil- 
ity between the two conjugate states. This geometrical 
picture gives an interpretation of the decoherence pro- 
cess as a measurement. The "measuring device" com- 
posed of nuclear spins evolves in two pathways, each for 
an electron spin eigenstate, up or down. When the sum 
of the distinguishability between the conjugate states 
of all pair-flips is large enough, the "device" states are 
"macroscopically" distinguishable. The coherence be- 
tween the basis states is destroyed by the measurement.!^ 
The pseudo-spin description will serve as the physical 
guide of the decoherence control by pulse sequences. 
III. FREE-INDUCTION DECAY 



The separation of the single-system dynamics from the 
ensemble dynamics offers an opportunity to study the 
FID of electron spin coherence due to dynamical en- 
tanglement, which would otherwise be concealed by the 
much stronger effect of inhomogeneous broadening. 

The pseudo-spin evolution is directly calculated from 



Eq. (29) to be 



kkit)) 



xtt 



icTk ■ sm 



I i), (33) 



where the sine function of a vector v is defined as 
sin(v) = (v/w) sin(i'). The distinguishability between 
the conjugate pseudo-spin states \tp^) is 



Si = {BkEk 



2 ,,:„_2Xa; ^■„2 



AkDk) sinc--|-smc 



^fc)cos^smc^ 



A \ Xh i ■ Xk ^ 
Ak) cos ^smc^ 



r 



(34) 



This expression is used in the numerical evaluation of the 
electron spin coherence. 

To gain some insight into the FID features, we make 
two approximations which are well justified. First, the 
pseudo-spins are separated into two groups, Qa and Gs, 
corresponding to non-local and local fiip-fiops, respec- 
tively. For the local pair-fiips {k G Gb), the transition 
matrix element is dominated by the intrinsic nuclear spin 
interaction as Bk ^ Ak under a moderate magnetic field, 
while for the non-local pair-fiips (fc e Ga)j the extrin- 
sic hyperfine-mediated interaction dominates. Thus the 



coherence is factorized as 



c 



A/I 



n li^'kim 




exp 



(35a) 



(35b) 



As shown in Fig. [6] the pseudo-spins in Ga and those in 
Gb have dramatically different precession behaviors. The 
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FIG. 6; (a) The conjugate Bloch vectors for a pair-flip by the 
intrinsic nuclear interaction with the diagonal part neglected 
[Xk ~ (Sfc, 0, ±i?fc)]. (b) The same as (a) but for a non-local 
pair-flip [x^ = ±{Ak,0, Ek)]. (c) and (d) are the projections 
of (a) and (b) to the x-y plane, respectively, in which the 
dotted trajectories show the effect of the diagonal nuclear 
interaction (Dfc) (with the deviation from the solid curves 
exaggerated for visibility). 



conjugate pseudo-spins for a non-local pair-flip precess 
into opposite directions, since both the transition ampli- 
tude (Afc) and the dominating part of energy cost {Ek) 
are associated with opposite signs for opposite electron 
spin states because of their hyperfine interaction origins. 
On the other hand, the intrinsic nuclear spin interaction 
is independent of the electron spin state, and the cor- 
responding pseudo-spin bifurcates tangentially. In the 
second approximation, we neglect the diagonal nuclear 
spin interaction Dk, which is justified from Eq. (34 1 with 
the condition Dk ^ Ek- 



This simplification can also be 
understood with the geometrical picture shown in Fig. |6] 
A small change of the energy cost by the diagonal nu- 
clear interaction induces only a slight modification of the 
pseudo-spin precession pathways, and hence a negligible 
change of the distance 5^. We will see, however, the 
diagonal nuclear spin interactions are important for the 
decoherence dynamics under control by pulse sequences. 

With the two simplifications and the conditions Ek ^ 
B],, A],, the distinguishability between conjugate pseudo- 
spin states is 



At^Al cos2 ^sinc^^^ « U-^Al [k e Ga) 

{k e Gb) 



t^ElBlsinc^^^t^ElBl 



, (36) 



where the short-time approximation holds for t <^ E^^ . 
Within timescales of interest {t <C B^^ ,A'^^), it is always 
satisfied that (5? <C 1, so the coherence is approximated 



_ exp I ~2t^ Al cos 
Cl_ - exp (-\t' ElBlsn.c'^\ 



2 ^sinc^^ , (37a) 



(37b) 



As the number of pseudo-spins is large, the details of in- 
dividual pseudo-spins are not important to the decoher- 
ence but what matters is the excitation spectrum, which 
is defined a^ 



keGA 

Ssie) ^ Y 5{e^Ek)Bl 



(38a) 
(38b) 



fcet 



for the pair-fiips in group Ga and G b , respectively. The 
decoherence in terms of the excitation spectra is 



>et 



'2 )' (39a) 



C^,^ = exp ^— y deSA{£) cos^ ^sinc 
Cl^ ^ cxp(^-^t* J deSB{e)e^smc^^^ . (39b) 



We are now ready to analyze some important features 
of the FID, including the short-time behavior, the depen- 
dence on the field strength and on the QD size, and the 
emergence of Markovian decay. 

In the timescale t <C B^^ , the pair-fiips cannot be 
described by energy-conserving scattering events, but 
should be understood in terms of quantum evolution. 
The decoherence is a highly non-Markovian dynamics. 
Particularly, in the very initial stage {t <C E^^), accord- 
ing to Eq. ( 36 1 , the electron spin coherence is well ap- 



proximated as 



'CI 



exp (-tVTf,^) , 

exp {-tyrls) , 



(40a) 
(40b) 



which is not an exponential decay, a typical indicator 
of Markovian dynamics. The decoherence times in the 
short-time limit are 




^1/2 .1/2- 



Here we have used the facts that Ak ^ N^^Aa^~^, 
the number of spin pairs connected by the extrinsic 
hyperfine-mediated interaction is about N'^, the number 
of spin pairs connected by the intrinsic interaction is in 
the order of N, and the typical hyperfine energy cost Ek 
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FIG. 7: (a) Electron spin coherence as lunctions of time for 
various field strengths, (b) The logarithm plot of (a), in which 
the curve for Bcxt = 12 Tesla is compared to the contribution 
by the hyperfine-mediated interaction (the dashed line) and 
that by the intrinsic interaction (the dotted line), respectively. 
The InAs dot is of the size 33 x 33 x 3 nm"^ and the nuclear- 
spin initial state \J) is randomly selected from an ensemble 
at temperature 1 K. The field strength is indicated by the 
numbers for each curve. 



FIG. 8: (a) The electron spin coherence in a QD of size 
15 X 15 X 2.4 nm^, under a field of 10 Tesla and tem- 
perature of 1 K. The lower curves are calculated with the 
hyperfine-mediated pair-flips neglected, and the dotted lines 
are the short-time profile. The non-Markovian-to-Markovian 
crossover is observed, (b) The excitation spectra for the non- 
local and local pair-fiips. 



for spin flip-flops between neighboring nuclei is in the 
order of f^-^l'^a^^^ - TV-^/a^^. 

A few features of the FID in the short-time limit 
(t <^ E^^) are summarized as follow. The decoherence 
caused by the intrinsic interaction has a quartic expo- 
nential decay profile and that by the hyperfine-mediated 
interaction has a quadratic one. So the decoherence is 
initially dominated by the hyperfine-mediated interac- 
tion and then cross over to the regime dominated by 
the intrinsic interaction as time increases. For a QD 
with 10^ nuclear spins and an external field of strength 
10 Tesla, the decoherence times due to the local and non- 
local pair-fiips, Tb and Ta, are estimated to be both of 
the order of 1 /is. The two decoherence times have very 
different dependence on the QD size (measured by N) 
and the field strength (measured by fig), as can be seen 
from Eq. (41 1. The hyperfine-mediate interaction could 



be the dominating decoherence mechanism when the QD 
size is small or when the external magnetic field is weak 
(e.g., Bcxt ~ 1 Tesla). The intrinsic interacton becomes 
dominating for large QD or when the field is strong (see 
Fig. [7]). InaQDof proper size and under an external field 
of proper strength, the crossover from the quadratic to 



the quartic exponential decay presents in a visible regime, 
i.e., in a time range where the electron spin decoher- 
ence has decayed by a finite amount but not vanished 
yet [see the curves for Bext = 12 Tesla in Fig.|7](b)]. The 
crossover can be tuned to occur after the coherence has 
vanished (as for Bext < 10 Tesla in Fig. |7| or before the 
decoherence is significant (as for B^Kt > 14 Tesla for the 
QD in Fig. [t] in which the data are not shown), then the 
decoherence in the visible regime will be dominated by 
the hyperfine-mediated interaction or the intrinsic inter- 
action, respectively. 

In general, the Markovian decay will emerge with time 
increasing, particularly for t > iJjT^, when the Fermi- 
Golden rule starts to come into effect [as indicated by 
the sine function in Eq. (37 1 and (39l]. Actually, the spin 



coherence at time t is mostly determined by the pseudo- 
spins with energy cost Ek < t~^, while the effects of 
pseudo-spins with higher precession frequency are can- 
celled out by destructive interference (between the fast- 
oscillating sine functions). In the long-time limit (while 
t ^ B^^ is still satisfied), the excitation spectra Sa/b{£) 
in Eq. ( 39 1 can be taken as "fiat" in the frequency range 
[-1/t, l/t], which covers those pseudo-spins that con- 
tribute significantly to the decoherence. The assumption 
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of "flat" spectra is just the Markovian approximation. 
So for time much greater than the inverse widths of the 
excitation spectra (which are roughly the inverse typi- 
cal energy cost of the pair-flips, E'^^), the electron spin 
coherence is approximated as 



£^A ^ ^—tSA{i / 4 cos^ a:sinc^a;dx 
£^B ^ / 4 sin^ a:sinc^a:(ia: 



e (42a) 
e"'^, (42b) 



presenting the exponential decay, a signature of the 
Markovian dynamics. The decoherence rate in the 
Markovian regime is determined by the weighted den- 
sity of states of the elementary excitations, i.e., the pair- 
flips, which, by the definition of the excitation spectra in 



Eq. ( 38 ) , is related to the Fermi-Golden rule by 



k 



(43a) 
(43b) 



where £,a/b denotes the integral in Eq. (42 1 and is of the 
order of unity. 

Under the conditions considered in this paper, the 
Markovian behavior will not be fully developed before the 
coherence has totally vanished, but instead a crossover 
behavior, namely, a reduction of the exponential index 
could be observed in the visible regime for certain QD size 
and field strength. An example of such non-Markovian- 
to-Markovian crossover is shown in Fig. [s] (a). The typ- 
ical energy cost of local pair-flips and that of non-local 
ones is in the order of AaN~^/^ and AaN^^, respec- 
tively, so the excitation spectrum of local pair-flips is 
much narrower than that of non-local ones [as shown in 
Fig. [s] (b)]. Thus we expect that the Markovian dy- 
namics will emerge earlier for the quadratic exponen- 
tial decay (caused by non-local pair-flips) than for the 
quartic one (caused by local pair-flips). To observe the 
non-Markovian-to-Markovian crossover before the deco- 
herence is complete, one should have E^^ < T2^a/b 
which is satisfied in relatively small QDs (for decoherence 
contributed by local pair-lips) or under relatively strong 
external field (for decoherence contributed by non-local 
pair-flips) . 

Due to the rich crossover behaviors discussed above, 
the decoherence in general cannot be characterized by a 
single time parameter T2. We introduce a decoherence 
time Tj^/g to quantify the time when the electron spin 
coherence has decayed to 1/e of its initial value. The 
decoherence time Ti/g can approach T2.A, T2,b, Too, a, or 
Too,B, depending on the QD size and the fleld strength 
which set the conditions for the crossovers. The field- 
and the size-dependence of different decoherence times 
shown in Fig. |9] and Fig. [TO] respectively, are consistent 
with the analysis above. For comparison, we also plot the 
decoherence time measured by the Hahn echo Tjf , which 
will be studied in the next Section. 
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FIG. 9: Field dependence of decoherence times (Ti/e - time 
for FID coherence being 1/e of its initial value, T2,a - FID de- 
coherence time resulting solely from hyperfine-mediated pair- 
flips, T2,B - FID decoherence time resulting solely from the 
intrinsic nuclear spin interaction, Th - decay time of the Hahn 
echo signal). The \/272,_b is plotted to compare with the Hahn 
echo decay time. The QD is as in Fig. [7] 
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FIG. 10: QD-size dependence of decoherence times (see Fig. [9] 
and the text for definition). The ^/2T2,b is plotted to compare 
with the Hahn echo decay time. The QD size is varied with 
fixed width:depth:height ratio 33:33:6, the field strength is 10 
Telsa, and the temperature is 1 K. 



IV. HAHN ECHO SIGNAL 

The FID studied in the previous section is not visi- 
ble in ensemble experiments in which dephasing due to 
inhomogeneous broadening is faster by orders of magni- 
tude than decoherence by dynamical entanglement. It 
might be difficult in the near future to directly observe 
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the single-system dynamics by filtering out the inhomo- 
geneous broadening w ith proiec tive measurement of the 
local Overhauser fieldA^ l ^^ l * Thus the spin echo may 
be used to study the electron spin decoherence problem 
in QDs. In Hahn echo experiments, the static Overhauser 
field experienced by the electron spin effectively changes 
its sign each time the electron spin is flipped by a short 
TT-pulse. Considering the simplest configuration in which 
only one pulse is applied at i = r, the precession phase 
accumulated from the random local fields is eliminated 
at < = 2t. Thus the decay of Hahn echo signals at t = 2t 
is solely due to the dynamical entanglement. The de- 
cay time of Hahn echo signals is usually used to quantify 
the spin decoherence time. We shall show, however, the 
Hahn echo decay time cannot be equated with the deco- 
herence time, as the decoherence by hyperfine-mediated 
interaction will be virtually eliminated from the echo sig- 
nal and that by intrinsic interaction is also suppressed to 
some extent. 

The TT-rotation or flip of the electron spin fo r Hahn 
echo can be operated by a GHz microwave pulse, '^^'^^'^ 
which is only marginally fast enough for eliminating the 
rapid dephasing by inhomogeneous broadening in III-V 
compound QDs. In optical control of an electron spin,^^ 
an arbitrary rotation of the electron spin can be com- 



pleted in the timescale of 10 ps via exciton- mediated 
Raman processes. With respect to the timescale of the 
electron spin decoherence, the optical pulses can be con- 
sidered instantaneous. The recent experiment on double 
GaAs QDs also employs a rather long DC voltage pulse 
to control the singlet-triplet transition to realize the spin 
echo,'^ which does not satisfy the instantaneous pulse 
condition. It is certainly interesting to study the elec- 
tron spin decoherence under the control of finite-duration 
pulses, but in this paper we would rather focus on the 
case of instantaneous pulses. 

As our general rule, the Hahn echo signal is deter- 
mined by the distinguishability between the conjugate 
pseudo-spins. In the rotating reference frame rested on 
the electron spin, the pseudo-spin states after the 7r-pulse 
applied at r are 

HUt>r))^U^{t-r)U^{r)\i), (44) 

which is equivalent to the conjugate pseudo-spin states 
IV'fc (t)) exchanging their pseudo-fields Xk when the elec- 
tron spin is fiipped. The trajectories of the Bloch vectors, 
projected to the x-y plane, are shown in Fig. |ll[ both for 
local and for non-local pair-fiips. The pseudo-spin states 
after the fiip pulse are 



COS cos ^ sm ■ sm 



cos ^ sin ^ 



Xt^ ■ Xp , . Xp . Xfr 

- cos — =^ sm — ^ h sm — ^ — x sm „ 



(45) 



The Hahn echo signal is 



(46) 



With the same arguments as in Sec. |III[ the electron 
spin decoherence can be factorized into two contribu- 
tions, namely £^ _ = £^ _ x _ with 



^A/B 



n expH^/2) 



(47) 



fceG 



A/B 



For the local pair-fiips driven by the intrinsic nuclear spin 
interaction, the distance at i = 2r is calculated by 
setting Dfc = Afc = 0: 



51 



]{2rfElBl 



(48) 



where the approximation in the second line holds for r <C 



Ey, ^ . The short-time behavior is 



£f^_(2T)«e-(2r)VT;^,B^ (49) 
the same as the FID signal except that the decay time 



H.B 



= V2r, 



2,B- 



(50) 



For non-local pair-fiips driven by the hyperfine-mediated 
interaction, the pseudo-fields are opposite to each 
other if the energy cost due to the diagonal nuclear spin 
interaction is neglected {Dk = 0). In this case, the 
pseudo-spins simply reverse their precession directions 
when the electron spin is fiipped, returning to its original 
states at the echo time, as shown in Fig. 11 (b), so the 



decoherence by the hyperfine-mediated pair-fiips is fully 
eliminated if is negligible. As the leading order of the 
distance 6k vanishes at the echo time, the second-order 
effect of the diagonal nuclear spin interaction becomes 
important. The distance 5k for non-local pair-flips, in- 
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eluding the efFeet of Dk , is given by 



corresponding to 



Si - \{2r)^DlAlsinc 
1 



EkT 



(51) 



where again, the approximation in the second hne holds 
for T ^ E^^- The residual decoherence due to the inter- 
play of the hyperfine-mediated interaction and the diag- 
onal nuclear spin interaction has the quartic exponential 
form for r ^ E^^ as 



The corresponding decoherence time is 



\keGA I 



1/2 



(52) 



(53) 



proportional to the square root of the field strength. 

Now we compare the echo signals to the FID signal 
in single-system dynamics (see Figs. [9] and 12). The ob- 
servation is that the decoherence due to the dynamical 
quantum entanglement in general is not measured by the 
Hahn echo signals. First, the FID and the echo signals 
could have qualitatively different short-time behaviours, 
as the former could have quadratic exponential decay due 
to the hyperfine-mediated interaction while the latter al- 
ways has the quartic exponential decay. Second, the de- 
cay timescale could be dramatically different from the 
FID to the Hahn echo. The difference in the decay time 
can be discussed in three regimes of field strength and 
QD size, roughly divided by 

I: ne<Ao.N-^^^, (54a) 
II : A^N-'/' « < Al/^N-y^^B-'/\ (54b) 

(54c) 



III : Al^'N-y^'B^'/' « n 
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FIG. 11: (a) The trajectories (projected to the x-y plane) of 
the conjugate pseudo-spins by the intrinsic interaction. The 
conjugate pseudo-spins l^'jl^) exchange their pseudo-fields Xk 
aX t = T when the electron spin is fiipped by a short pulse. 
The solid and dotted curves denote the trajectories before and 
after the flipping pulse, respectively, (b) The same as (a) but 
for a non-local pair-flip by the hyperflne-mediated interaction. 



II: Th,a:S>Th^b 
III: r2 4>r2R. 



T2,B > T, 



2.A, 



(55a) 
(55b) 
(55c) 



In regime I, the hyperfine-mediated interaction is im- 
portant both in the FID and in the Hahn echo decay. 
As the excitation of non-local pair-fiips is eliminated in 
the leading order at the echo time [See Fig. 11 (b)], the 
echo decay time is much longer than the FID decoherence 
time, namely. 



T2,a/Th,a - VMn^Dk/Al « 1, 



(56) 



for QDs of reasonable sizes. The pair-correlation ap- 
proximation, however, is not valid in regime I, since the 
condition that the number of pair-flips in the timescale 
of Th,a be much fewer than the number of nuclei re- 
quires N ^ Ai/^fle^^'^D, ^^"^ for the hyperfine mediated 



non-local pair-flips and N <C Ai/'^fle^^^D'^^^'^ for the lo- 
cal pair-flips driven by the intrinsic interaction (see Ap- 
pendix |C] for estimation of the pair-flip numbers) , which 
cannot be simultaneously satisfied. Theories beyond the 
pair-correlation approximation need to be developed to 



explore this regime. In regime II defined by Eq. (54b I, the 



hyperfine-mediated interaction contributes significantly 
to the FID but virtually nothing to the Hahn echo decay. 



FID 

FID, local pair 
Hahn echo 




FID 

FID, local pair 
Hahn echo 



2t ([Xs) 



FIG. 12: (a) Comparison of the Hahn echo (dashed green 
line) and the FID (solid red line) signals. The FID signal is 
also shown with the hyperfine-mediated pair-fiips neglected 
(dotted blue line), (b) The logarithm plot of (a). The QD is 
as in Fig. [7| with Bext = 2 Tesla. 
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In this regime, the decay time measured in echo signals 
is much longer than the FID decoherence time as 



T, 



= V2T2,B > T, 



(57) 



In regime III, the hyperfine-mediated interaction is 
strongly suppressed, so the FID and the Hahn echo decay, 
both determined by the intrinsic nuclear spin interaction, 
have essentially the same decay profile but different de- 
cay times (differing by a factor of v^). Thus, in all the 
three regimes, there is substantial difference between the 
FID and the Hahn echo decay. Such difference is in- 
deed a consequence of the modified dynamics of a meso- 
scopic bath (the nuclear spins) under the manipulation of 
the quantum object (the electron spin). In semiclassical 
spectral diffusion theories, the electron spin experiences 
passively the "background" of fluctuating local fields. In 
the full quantum theory, the electron spin dynamics ac- 
tively alters the mesoscopic bath dynamics (as the nu- 
clear spin pair-flips depends dramatically on the electron 
spin state). Such a view paves the path toward manip- 
ulation of the nuclear spin dynamics in mesoscopic QDs 
and control of the electron spin decoherence. 

The active modiflcation of the bath dynamics by the 
manipulation of the quantum object also manifests it- 
self in the effect of the diagonal nuclear spin interaction 
on the Hahn echo signal. The diagonal nuclear spin in- 
teraction contributes much less energy cost of a pair-flip 
than the hyperfine interaction (Dk <^ Ek) and has no 
direct effect on the effective Overhauser field felt by the 
electron spin. Thus it is expected that the diagonal inter- 
action terms should have little effect on the the electron 
spin decoherence. The calculation of the FID signal in 



Sec. Ill indeed confirms such a semiclassical argument. 
The difference between the quantum theory and its semi- 
classical counterparts becomes significant when the elec- 
tron spin is under control. In semiclassical theories, the 
electron spin only passively detects the fiuctuating local 
field which is essentially independent of the electron dy- 
namics, so the diagonal nuclear spin interaction should 
contribute no more than it does in the FID configura- 
tion. In the quantum theory, the nuclear spin dynamics 
conditioned on the electron spin state is held responsible 
for the electron spin decoherence. The diagonal nuclear 
spin interaction affects the nuclear spin dynamics signifi- 
cantly and thus contributes to the decoherence when the 
nuclear spin dynamics is under active control by the ma- 
nipulation of the electron spin state. In next sections, 
we will see the importance of the diagonal nuclear spin 
interaction in determining some qualitative decoherence 
features. 



V. DISENTANGLEMENT AND 
RECOHERENCE 

In the previous section we have seen that the decoher- 
ence caused by the hyperfine-mediated interaction can 



be nearly eliminated in Hahn echo signals. Such elimi- 
nation of decoherence is realized because the hyperfine- 
mediated pair-fiips are reversed when the electron spin 
is flipped aX t — T [see Fig. 11 (b)] and hence the 



pseudo-spin states |'0fc^(^)) return to their original po- 
sitions &i t ~ 2t. As the decoherence is caused by 
the entanglement between the electron spin and the nu- 
clear spins, the elimination of decoherence, also called 
"recoherence" , can be understood in terms of disentan- 
glement. For non-local pair-flips under Hahn echo con- 
trol, the disentanglement is realized at t = 2t when the 
electron spin and the pseudo-spin are brought back to 
a factorized state C+\+)\ip^ {2t)) + C^\-)\ip^ {2t)) with 
|V'±(2r)) = |V±(0)) = U). 

In general, a pseudo-spin (or the environment) is dis- 
entangled from the electron spin (or the quantum sys- 
tem in contact with the environment) whenever its states 
\ip^{t)) for different electron spin states |±) become iden- 
tical, but need not return to the original states |V'fc^(0))- 
When the electron spin is flipped, the conjugate pseudo- 
spins exchange their pseudo-fields, so their trajectories on 
the Bloch sphere will inevitably intersect, leading to the 
disentanglement.'^ The intersection of conjugate pseudo- 
spins driven by the intrinsic nuclear spin interaction is ac- 
tually seen in Fig. 11 (a), which, unlike that for non-local 
pair-fiips, occurs at a time different from the echo time. 
The intersection time, in general, is different for differ- 
ent pair-flips, making it impossible to fully recover the 
lost electron spin coherence. Nonetheless, in the short- 
time limit (r <C E^^), the distance between conjugate 
pseudo-spins can be eliminated in the leading order of 
time, simultaneously for the same group of pair-flips, as 
have been seen for non-local pair-flips. It is straightfor- 
ward to verify that for the local pair-flips, such coinci- 
dent disentanglement does take place at t = \f2T in the 
single-pulse Hahn echo conflguration [the magic number 
\p2 can be understood by noticing that oc at short 



With the leading order contribution vanishing, the dis- 
tinguishability between conjugate pseudo-spins at the 
disentanglement time (\/2t), including the effect of the 
diagonal nuclear spin interaction, in the next leading or- 
der of T for local pair-flips is 



4 (V27 



EkBkDkT^ 



(58) 



Note that the residual decoherence would otherwise be 
of a higher order in r if Z?j. were set to zero: 



(59) 



Here we see again that the small effect of the diagonal 
nuclear spin interaction (the D^. term) emerges when the 
nuclear spin dynamics is under control. 

A semiclassical spectral diffusion theory also predicts 
a coherence recovery at a time earlier than the echo time 
(2t) when the random field memory time (tj,) is compara- 
ble to the pulse delay time.^^ But the semiclassical theory 
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differs qualitatively from the quantum counterpart in at 
least two aspects: (1) The recovery time depends on t^, 
being about r + for < t and approaching 2r when 
Tc ^ t; and (2) The exponential index of the decoher- 
ence profile in the spectral diffusion theory is essentially 
unchanged by the control. The latter difference points 
directly to a fundamental shortness of the semiclassical 
picture: It does not take into account the active control 
of the local field fluctuation by the flip of the quantum 
object but considers the quantum evolution of the bath 
as a fluctuating field experienced passively by the quan- 
tum object. 

Inclusion of contributions from both local and non- 
local pair-fiips fields the pulse-controlled decoherence in 
the short-time limit, 

where the effective decoherence time due to the residual 
entanglement from local pair-flips is deflned by 



6 ^ V 

off,B / . 



(61) 



The electron coherence is recovered at V2t, even when 
the pulse is applied after the coherence has been totally 
lost in FID (for r > T2.b)- Such recoherence in single- 
system dynamics and the spin echo in ensemble dynam- 
ics have different physical bases: the former is due to 
quantum disentanglement, but the latter is due to clas- 
sical refocusing of random phases. The difference be- 
tween the two is already evidenced by their different oc- 
currence time. Fig. 13 (a) plots an example of the real- 



time dependence of the electron spin coherence under a 
single-pulse control, which demonstrates the recovery of 
the coherence even when the FID signal has vanished 
for the chosen delay time. The contribution from the 
hyperfine-mediated pair-flips, which are not reversed at 
•\/2r, makes the coherence to be only partially recovered. 
Under a stronger magnetic field, a full recovery at \/2r 
can be realized as the hyperfine-mediated interaction is 
fully suppressed (not shown). 

(a), 



As illustrated in Fig. 
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a sequence of pulses can 
force the trajectories of conjugate pseudo-spins to cross 
into each other again and again. In the short-time limit 
(the delay time r <C E'j^^), an equally spaced pulse se- 
quence eliminates the entanglement due to the local pair- 
fiips up to the terms, leading to a sequence of reco- 
herence at 



1 r, 



1,2, 



(62) 



One can verify these magic numbers by using the 
quadratic dependence of Sk on time and checking that 

(^^yn{n+l)T^ - (nr)2 = (nr)^ - (^^(n- l)nr^ . The 

decoherence caused by the hyperfine-mediated pair-fiips 
under the control of evenly spaced pulses is bounded, 
which can be understood from the reversed precession of 
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FIG. 13: (a) The electron spin coherence under the control 
of a short 7r-pulse applied at r = 2 /^s (when the FID signal 
has vanished), the recoherence at \/2t is pronounced while 
no signal survives at the echo-time 2r. (b) The electron spin 
coherence under the control of a sequence of evenly spaced 
pulses, (c) The electron spin coherence under the control of 
a Carr-Purcell pulse sequence. The arrows indicate positions 
of the short pulses. The solid blue (dotted red) lines are 
calculated with(out) including the hyperfine-mediated pair- 
fiips. The QD is as in Fig. [7| with Boxt = 10 Tesla. No 
ensemble average is done. 



the corresponding pseudo-spins. Fig. 13 (b) confirms our 
above analysis. By using a long sequence of pulses at 
short intervals, the electron spin coherence can be pre- 
served for an arbitrarily long time until other mechanisms 
of decoherence (such as phonon scattering) come into ef- 
fect. 

The recoherence in the control configurations discussed 
above, however, is not observable in ensemble experi- 
ments, since the rapid dephasing due to inhomogeneous 
broadening will prevent any coherence being observed ex- 
cept around the spin-echo time. It is desirable to design a 
pulse sequence to force the recoherence from disentangle- 
ment to coincide with a spin-echo from phase-refocusing 
so that the disentanglement effect can be studied in en- 
semble experiments. This possibility can be seen from 
Fig. 14 (b). When the electron spin is flipped by two 



short pulses, for example, the disentanglement time after 
the second pulse can be adjusted by changing the delay 
time between the two pulses. It is shown straightfor- 
wardly that when the pulse sequence is designed such 
that the first pulse is applied &t t = t and the sec- 
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FIG. 14: Trajectories of conjugate pseudo-spins from the in- 
trinsic nuclear spin interaction (a) under the control of a se- 
quence of pulses equally spaced and (b) under the control of a 
Carr-Purcell two-pulse sequence. The green bars on the time 
axis represent the short pulses flipping the electron spin. 



ond at t = 3t, the disentanglement coincides with the 
echo time at t ^ At (by virtue of the change in Sk, 



(4r)2 - (3r)2 - [(3r)= 



0). Obviously, the 



hyperfine-mediated non-local pair-flips are disentangled 
from the electron spin at the echo time. A two-pulse 
sequence so designed happens to be the famous Carr- 
Purcell sequence, widely used in NMR experiments to 
dynamically decouple the nuclear spins. The disentan- 
glement studied here, however, is fundamentally different 
from the dynamical decoupling, as will be discussed later 
in the next section. 

In the rotating frame of the electron spin, the pseudo- 
spin states at the echo time after a Carr-Purcell sequence 



|V,^(4r)) 



i) 



i) 



r/2| 



i), 



with the series of angles of rotation defined by 



t/2 



(63) 

(64a) 
(64b) 



The distinguishability between conjugate pseudo-spins 
at the echo time At in the leading order of delay time is 

n2 



51 



(66) 



for T <^ Ef. . Note the decoherence would otherwise be 



of a higher order in r if is set to zero, 
51 - 2^\^BlEl. 



(67) 



This demonstrates again the role of the diagonal nuclear 
spin interaction. We note that the decoherence profile in 
the two-pulse control configuration, in particular the ex- 
ponential index, is qualitatively different from the predic- 
tions of semiclassical theories^ in which the bath dynam- 
ics is accounted by a randomly fluctuating force (without 
being actively changed by the manipulation of the quan- 
tum object) and the resultant decoherence exponential 
index is basically unchanged by the pulse control. The 
recoherence at the echo time 4t in the Carr-Purcell con- 
trol configuration is clearly demonstrated in Fig. [13] (c), 
even though the signal at the Hahn echo time 2t after 
the first pulse is absent. 



VI. CONCATENATED DISENTANGLEMENT 

A. Formalism and geometrical picture 

We have used the iteration of effective precession an- 
gles to derive the pseudo-spin states under the Carr- 
Purcell control in Eqs. ( 63p5 1. Such an iteration for- 
malism inspires us to borrow the idea of concatenation 



control in dynamical decoupling recently studied for jjre- 
serving coherence in quantum computation. ^-'^^^ To 
illustrate the basic idea, we can imagine that if Xk 



Eq. (64a I is replaced by the effective field corresponding 



to the precession angles at the Carr-Purcell echo time, 
i.e., ©J/t, the distance between the conjugate pseudo- 
spins would be eliminated to an even higher order of the 
delay time. Such a replacement can be done recursively, 
and a concatenated pulse sequence can be so designed to 
control the decoherence. The unit cell of pulses to be con- 
catenated does not have to be the Carr-Purcell sequence, 
but it could be either simpler or more complicated, or 
even various types of unit sequences can be interwoven 
to construct a sophisticated concatenation. 

To exemplify the concatenated disentanglement for de- 
coherence control, we choose the simplest concatenation 



The vector Of has a geometrical interpretation as the 
pseudo-spin effective precession angle 9'^ about the axis 

along ef at t, 2t, and 4r for Z = 0, 1 and 2, respectively sequence which is recursively defined by Eq. ([64|. Fig.[l5 
(the pseudo-spin index k is understood where no confu- 
sion is caused) . Eq. ( 64b I yields a geometrical recursion 
relation: 



e 



sm 



cos 



3^ 



■ -I 

cos sm — ^ 
2 2 



9i ■ Of 

cos sm — ^ 

2 2 



e 



■ Of . ^, 

' sm — ^ X sm — 
2 2 



± 



cos COS — 

2 2 



± 



ef 



(65a) 
(65b) 



illustrates how such a sequence can be concatenated to 
any order. The 0th order is a FID evolution with no pulse 
control and the (Z + l)th order sequence is constructed by 
two subsequent /th sequences with or without one extra 
fiipping pulse inserted in between depending on whether 
I is even or odd. Accordingly, the evolution propagator 
is iterated as 



^-txi-^kT/2 ^ -i0^-&/2 



•<t/2 



(68a) 
(68b) 
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FIG. 15: Concatenated sequences of short-pulses flipping the 
electron spin, represented by vertical bars. 



For a pictorial understanding of the concatenated con- 
trol of decoherence, we rewrite the effective precession 
angles as 







Ri±r 



(69) 



Without confusion, the vectors will also be referred to 
as effective precession angles since « 0^/2 for small 



precession angles. As depicted in Fig. 16 (a), the conju- 
gate precession angles are decomposed into the common 
part R; and the difference part r;. Then, the recursion 



in Eq. ( 65 1 is rewritten as 



To 
Ro 
Ri 

R( 



1 . xIt 1 . XfcT 

- sm — sm , 

2 2 2 2 ' 

1 • Xfe-r , 1 . XfeT 

- sm — — sm — — , 

2 2 2 2 

COS — — sm — h cos — ^ sm — — , 

2 2 2 2' 

2R,_ixri_i, (?>1), 



— 2R;_i w 1 — Ri_i 



a>2). 



(70a) 

(70b) 

(70c) 
(70d) 
(70e) 



The distance 5^ is in the same order of the difference 
part. Particularly, for the /th order concatenation, the 
distance at r; = 2't is 



51 



1 — rij + • cr 



l-nf 



= 1 - 



Art 



1 — ?if + rij • — i2ri ■ a J 1 — — i (n^ x n 



l-rf- \zi+iRi-ziJ\~Ri 



(71) 



for / > 1, where ri; — — — ~\- rf and Zi is the 

z-component of Vi/ri. 



B. Short-time behaviors 

To gain some insight into the concatenated control of 
decoherence, first consider the short-time limit r <C . 
There, the recursion of the rotation angles has an intu- 
itive form as 

Ri « 2'Ro^2V ^^ + ^" ^(2i?fc,0,jJfc)|, (72a) 

ro « (x^-Xfe)^ = (2Afe,0,i?fc)^, (72b) 
« 2'RoXri_i, (Z > 1), (72c) 

until the precession angle Ri approaches the order of 
unity at a threshold order of concatenation /q given by 

/o«-log2(Bfcr). (73) 



The distinguishability between conjugate pseudo-spins in 
the leading order of delay time is 

<5f - 4rf (1 - ) , (74) 

determined by the difference between the conjugate pre- 
cession vectors. The concatenated control can be under- 
stood in the geometrical picture shown in Fig. [l6] (b): 
In the FID (the 0th order concatenation), the common 
part of the precession angles i?o is roughly determined 
by the nuclear spin interaction strength times the time, 
i.e., B^T, which is much less than 1 in the timescale of 
interest, and the difference part ro is roughly the pre- 
cession angle given by the hyperfine energy cost E^. In 
the 1st order concatenation, the common part Ri is ap- 
proximately in the same direction of Rq with the angle 
roughly doubled if E^t ^ 1, and the difference part ri 
is perpendicular to both Rq and Tq with the amplitude 
reduced by a factor of 2i?o from tq. Staring from the 
2nd order concatenation, the common part R/ will be 
along the the same direction as the 1st order one Ri and 
the difference part r/ will be alternatively in the two or- 
thogonal directions ri and V2, both perpendicular to the 
common part direction R;. By each further level of con- 
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catenation, the common part is increased by a factor of 
two and the difference is reduced by a factor of 2Ri for 
I < Iq. Thus each level of concatenation will reduce the 
difference between effective conjugate precession angles 
by an order of delay time times the nuclear spin inter- 
action strength and remove the decoherence accordingly, 
until the controlling effect is saturated at the threshold 
level Iq. 

With the short-time approximation of the precession 



angles given in Eq. ( 72 1 , the distinguishability between 



conjugate pseudo-spins in the leading order of r is 



{ABl + Diy \BuEk-AkDk)\ (75) 



with 



(/==1,3, ...) 

ABl{ABl + Dl)-' (/ = 2,4,...) ' ^^^^ 



Alternatively, for a fixed time t, a certain order of con- 
catenation control with t = gives 



51 



4^2/+22-i('+3) H 



{ABl + Diy ' (BkEk - AkDkf, (77) 



which can be suppressed to an arbitrary power of the 
echo time t times the nuclear spin interaction strength 
under the condition 

Ri ^ Bkt « 1. (78) 
The electron spin coherence under control is 



(79) 



decaying exponentially with powers of r at 2, 4, 6, 8, . . ., 
for the concatenation order 0, 1, 2, 3, . . ., respectively. 

To show the effect of the diagonal nuclear spin interac- 
tion in the decoherence control, we note that in Eq. ( 76 ) 




2:2m = 1 if Dk = and thus Sj. vanishes in the given order. 
In this case, the distinguishability between the conjugate 
pseudo-spins in the next leading order for I = 2m is 



51 



-42i(i+3)^2^2i+2 



(80) 



So when is put to zero, the coherence decays expo- 
nentially with powers of t at 2, 4, 8, 8, 12, 12, ... for 
concatenation order 0, 1, 2, 3, 4, 5, . . ., respectively. 

The effective preservation of the coherence by concate- 
nated control is demonstrated by the numerical simula- 



tions presented in Fig. 17 Note that the coherence shown 
in the figure is at the time 2't which doubles at every 
order of concatenation. The role of the diagonal part 
of the nuclear spin interaction is demonstrated [compare 



Fig. 17 (b) and (c)]. 



C. Controlling small parameters 



The suppression of the decoherence by concatenated 
control, however, does not depend on the short-time con - 
This point can be seen from Eq. ( 70 1 : 



dition T <^ Ef.^ 




T (^IS) 



FIG. 16: (a) The effective precession vectors for conjugate 
pseudo-spins are decomposed into the common component 
and the difference part, (b) The iteration of the effective 
precession vectors from the Zth order to the {I + l)th order 
concatenation, (c) The pseudo-fields are decomposed into the 
perpendicular and the parallel components. 



FIG. 17: (a) The electron spin coherence under the Ith order 
concatenation control, I = 0, ... ,5, as functions of the pulse 
delay time r. (b) The logarithm plot of (a), (c) The same 
as (b), but with the diagonal nuclear spin interaction put to 
zero. The QD is as in Fig. [13} 
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An iteration of concatenation suppresses the decoher- 
ence further as long as the common component of the 
conjugate precession angles i?/ ^ 1, without requiring 
EkT <^ 1. Since Rq is in the order of nuclear spin inter- 
action strength times the delay time (Bfcr), the condition 
for a concatenated sequence to be efficient in eliminating 
the decoherence is that the delay time be much shorter 
than the inverse strength of the nuclear spin interaction, 



r<^B7:\D-\A, 



(81) 



which is a much less stringent condition than t ^ i^jT . 
To calculate the electron spin coherence under this re- 
laxed "short-time" condition, we decompose the conju- 



gate pseudo-fields in such a way that 
Xfc - a - c. 



(82a) 
(82b) 



with a _L b II c, as depicted in Fig. 16 (c). Then, we 
are able to separate different timescales according to the 
orders of magnitude of different components. We notice 
that a and |6 — c| are of the order of nuclear spin inter- 
action strength and b and c are of the order of hyperfine 
energy cost Ek- Within timescales given by Eq. (81 1, the 



precession angles in the leading order of the nuclear spin 
interaction strength are 



To = sm 



n 



(h + c)T 



„ (b — c) T (6 + c) r ar (&+c)t 
Ro = ^ ' cos ^ +— smc ^ ^ , 

= 2'-2(b-c)r + 2'-iarsinc^^^^^^, (/ > 0), 



2'-^r'+^\EkBk^AkDk\{4Bl + Dl)'^ , {l>0), 



4 

i-i 



2 • 2 

a smc 



(6-f c)t 



(83a) 
(83b) 
(83c) 

(83d) 



(83e) 



under the condition Ri 
direction of r; is 

f 



<C 1. The z-componcnt of the the threshold concatenation level is estimated as 



4agsinc^(£fcT) 



2m - 
2m 



1) 
2) 



(84) 



The distinguishability between conjugate pseudo-spins 
is determined by Eq. (71 ). The decoherence is estimated 
accordingly. 

It is obvious that 
the maximum delay 
trol to be efficient 



strength ^/4El 



Dl 



the small parameter controlling 
time for a concatenation con- 
is the nuclear spin interaction 
instead of the object-bath inter- 
action strength (i.e., the hyperfine energy cost Ek). The 
threshold concatenation level Iq, beyond which an addi- 
tional level of concatenation cannot suppress the deco- 
herence further, is determined by the condition that the 
precession angle approaches 1, i.e., Ri^^ < 1. With the 
precession angle given by 



Ri 



2'^ V J4B2sinc^ {EkT) + D 



< 2^-^r,lABl + Dl 



(85) 



'l0g2 



4B2 



Dl 



(86) 



also determined by the nuclear spin interaction strength. 
A close examination of Fig. 17 (a) reveals that the deco- 



herence suppression is still effective even when the pulse 
delay time r is substantially longer than the typical value 
of the inverse hyperfine energy cost E'^^ , which is about 
0.5 ^s for the QD considered. 



D. Coherence stabilization 

The concatenated pulse sequence can be repeated. The 
propagator of a pseudo-spin for n repeated Ith-order con- 
catenation sequences is 



Ln 



{ufy = exp {-inOf ■ a/2) 



(87) 



Note that in implementing the pulse sequence for re- 
peated concatenation, a 7r-pulse should be inserted be- 
tween two unit sequences if / is an odd number. The 
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common and difference parts of the conjugate precession 
angles R/^„ and Ti^n can be defined similarly to those 
for unit concatenation sequences. The precession vectors 



9 



± _ 

l,n 



lOi do not changes their directions with increas- 



ing the number of repeated units but only the angle of 
rotations grows linearly with n. Therefore, 



ri,n/Ri,n = ri/Ri = 



Since Ri 



/2 < 1, the distinguisha- 



bility between conjugate pseudo-spins given by Eq. (71) 
is bounded from above by 



(89) 



Since by each order of concatenation until the thresh- 
old level is reached, the common precession angle Ri is 
increased by a factor two and the difference part r/ is 
suppressed by a factor of i?/ ^ 1, the ratio ri/Ri can 
be made arbitrarily small. In the mesoscopic system, 
the number of excitation modes (pseudo-spins) is finite, 
so the coherence under repeated concatenation control is 
bounded from below by 



C+- {uTi) > exp —2 



/Rf - £+,_ (rj_i) . (90) 



The summation 2rf / Rf can be made arbitrarily small 
by choosing a small delay time and a sufficiently deep 
concatenation level, which means the electron spin co- 
herence is stabilized virtually without decoherence. The 
coherence after n units of Ith order concatenation se- 
quences is no less that that after a single unit of (/ — l)th 
order concatenation sequence. The analysis above is ver- 



ified by numerical simulations presented in Fig. 18 
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FIG. 18: The electron spin coherence under the control of 
repeated units of concatenated sequences, as functions of the 
number of units for a fixed pulse delay time r = 0.2 fis. The 
QD is under the same condition as in Fig. |17| 



VII. SUMMARY AND DISCUSSIONS 

We have studied the electron spin decoherence in a 
mesoscopic bath of nuclear spins in semiconductor QDs. 
For QDs of interest in experiments, the electron- nuclear 
hyperfine interaction is much stronger than the nuclear 
spin interactions, so the electron spin and those nuclear 
spins in direct contact with the electron can be well iso- 
lated from the environment as a closed quantum system 
in the timescale of interest, giving the definition of a 
mesoscopic system. 

When the longitudinal spin relaxation is absent (sup- 
pressed here by a moderate to strong magnetic field 
^ 0.1 T and a low temperature ^ 1 K), pure dephas- 
ing (or transverse relaxation) induced by the nuclear spin 
bath is the dominant cause for the lost of electron quan- 
tum coherence. In the quantum theory, decoherence is 
a consequence of the entanglement between the quan- 
tum object and the bath established by the evolution of 
bath along bifurcated pathways in the Hilbert space for 
different basis states of the quantum object. By con- 
trast, semiclassical spectral diffusion theories ascribe the 
decoherence to accumulation of a random phase from a 
fluctuating local field. The nuclear spin dynamics, which 
is conditioned on the electron spin basis states, can be 
shepherded by manipulation of the electron spin. In par- 
ticular, when the electron spin is flipped, the bifurcating 
nuclear spin pathways will exchange their evolution di- 
rections in the Hilbert space and could intersect at a later 
time, leading to the disentanglement of the electron from 
the nuclei and hence the restoration of the coherence even 
after it has been totally lost to the nuclear bath. Such re- 
coherence by disentanglement is fundamentally different 
from spin echoes which result from re- focusing of random 
phases. 

The solution to the many-body nuclear bath dynamics 
relies on four essential conditions: (1) The QD size is in 
the mesoscopic regime. (2) The external fleld is much 
stronger than the hyperflne interaction. (3) The temper- 
ature is low for the electron spin but high for the nuclear 
spins. (4) The timescale considered is much shorter than 
the inverse nuclear interaction strength. Such conditions 
are realistic. They make possible two important simpli- 
fications in theory. First, the electron spin dephasing in 
ensemble dynamics can be factorized into two factors, 
namely, one due to the inhomogeneous broadening, and 
the other due to the entanglement between the electron 
spin and the nuclear spins. Second, the nuclear spin dy- 
namics in the timescale of interest can be approximated 
as independent excitations of pair-wise flip-flops. The so- 
lution of the bath quantum dynamics provides a simple 
geometrical basis for the design of pulse sequences for 
recoherence. The controllability of the mesoscopic bath 
dynamics, however, is a rather general concept indepen- 
dent of the simplifications made in the present paper. 

A close examination of the single-system dynamics free 
of the inhomogeneous broadening reveals a wealth of phe- 
nomena in the electron spin decoherence. For example. 
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a coherence recovery at a time different from the spin 
echo time would have been totally eclipsed in ensem- 
ble dynamics, since the dephasing due to inhomogeneous 
broadening is usually faster by orders of magnitude than 
the decoherence due to the entanglement. It is also shown 
that the single-pulse Hahn echo signals do not measure 
the electron spin decoherence in FID configuration: First, 
the decoherence caused by hyperfine-mediated nuclear 
pair-flips has a strong effect on the FID when the ex- 
ternal field is not too large, but it is removed from the 
Hahn echo signal due to disentanglement; Second, the 
decoherence by the intrinsic nuclear spin interaction is 
also partially suppressed, which results in a Hahn-echo 
decay time approximately \/2 times the FID time due 
to the intrinsic interaction. The decoherence behavior 
in single-system dynamics (such as the recoherence at 
•\/2t) could eventually be directly observed in experi- 
ments if the inhomogeneous broadening can be filtered 
out, e.g., by a project ive measurement of the local Over- 
hauser ^eldm^SmM 

To observe the recoherence by disentanglement in en- 
semble experiments,'^*^ pulse sequences can be designed 
to force the disentanglement to coincide with a spin echo. 
The simplest solution is the Carr-Purcell sequence, which 
has one pulse at r and another at 3t and forces the dis- 
entanglement to take place at the echo time 4r. More 
sophisticated designs such as concatenated sequences can 
be employed to suppress the decoherence to an arbi- 
trary order of the pulse delay time in units of the in- 
verse nuclear interaction strength, and a repetition of 
the concatenated sequences can preserve the electron spin 
coherence until other decoherence mechanisms (such as 
phonon scattering) become significant. 

The design of concatenated pulse sequences for 
disentanglement resembles the dynamical decoupling 
schemes which have been developed in decades for NMR 



spectroscopics 
computation 



45 47 48 
51 52 54, 



and recently studied for quantum 
The disentanglement schemes ex- 
ploit the controlling power of pulse sequences in a subtly 
different way from the dynamical decoupling schemes: 
The latter uses short pulses to flip the quantum object 
frequently so that the object-bath interaction is dynami- 
cally averaged to zero, while the former does not seek to 
eliminate the effective coupling but rather focuses on con- 
trolling the evolution of the bath wavefunction. For ex- 
ample, for local nuclear spin pair-flips in the single-pulse 
control conflguration, the effective interaction (8) HcS 
at the disentanglement time V^t, defined by HcS = 
H+g ~ and e-'"^fi^^ = e-'^^i^-^^^-e''"^^, does 
not vanish even in the leading order of the hyperfine cou- 
pling. This is also the case for the disentanglement under 
the control of equally spaced pulses. One may question 
why a non-zero effective object-bath interaction could 
lead to vanishing decoherence. The reason is that the 
control steers the wavefunction evolution and the initial 
state of the bath is of some special form (namely, the nu- 
clear spin states at the temperature under consideration 
possesses no off-diagonal coherence). The disentangle- 



ment time is in general dependent on the initial state 
(a universal disentanglement scheme independent of the 
initial state is only possible when the effective system- 
bath coupling vanishes). The dependence of the disen- 
tanglement on the initial state is implied in the control 
by evenly separated pulses shown in Fig. 14 (a): For 



instance, if the nuclear spins there at \/2t is set as the 
initial state (which possesses off-diagonal two-spin coher- 
ence), the disentanglement after the electron spin flip at 
To would occur at (\/3 — 1)(\/2+1)t'o (on the contrary, for 
initial states J') which has no off-diagonal coherence, the 
disentanglement occurs at \/2tq). It would be very inter- 
esting to study the decoherence and the disentanglement 
for spin baths with certain initial off-diagonal coherence 
such as in quantum memory applications.'^ 

Even though the concatenated disentanglement does 
cause the effective electron-nuclear coupling vanish, the 
difference from the dynamical decoupling is still unam- 
biguously evidenced by their different controlling small 
parameters. The controlling small parameter in dynam- 
ical decoupling is the object-bath interaction strength, 
and the pulse delay time should be much shorter than 
the inverse of that strength. In the dynamical disentan- 
glement, the object-bath interaction strength is a quite ir- 
relevant parameter, but the short-time condition is given 
by the bath interaction. Furthermore, in concatenated 
dynamical decoupling, the decoherence is eliminated in 
orders of the system-interaction strength times the pulse 
delay time, but in the disentanglement, the decoherence 
is suppressed in orders of the bath interaction strength 
times the delay time. 
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APPENDIX A: FACTORIZATION OF 
SINGLE-SYSTEM DYNAMICS AND 
INHOMOGENEOUS BROADENING 

The nuclear spin dynamics is basically determined by 
the excitation spectra of the pseudo-spins when the bath 
size is large enough. The excitation spectra for differ- 
ent nuclear spin configurations \J') should be the same 
except for a relative variance in the order of 1/ \/N ac- 
cording to the central limit theorem in statistics. Here we 
give a numerical test of the single-system dynamics for a 
few initial nuclear spin configurations randomly selected 
from the thermal ensemble. The results are shown in 
Fig. [19] Remarkably, the coherence for four initial states 
randomly selected from the thermal ensemble is visually 
indistinguishable, both in the FID regime {t < 0.5 /is) 
and after the controlling pulses are applied {t > 0.5 fis). 
The relative deviation between different curves is consis- 
tent with the estimate by the central limit theorem. 
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FIG. 19: (a) The electron spin coherence under the control 
of a Carr-Purcell pulse sequence, for various initial nuclear 
spin states (indicated by arrows in Fig. |4|: |j7i) with local 
Overhauser field Ej-^ = — 0.5r2 (dot-dashed red line), {^2) 
with Ej^ = (sohd black line), jj'a) with Ej^ = CSFj (dot- 
ted blue line), and ijTi) with Ej^ = F2 (dashed green line). 
The four lines are indeed indistinguishable in the figure, (b) 
The deviation of the coherence for various initial states from 
that for the state A£+,_(f) = \{Jn\e'"-'e-'"+'\J„)\ - 

\{Ji\e'"-'e- ■ 



\{Jn\e'"-'e-'"+'\ 
\Ji)\, ampHfied by 1000 times. The InAs 



QD is of size 34 x 34 x 3 nm'^, at a temperature of 1 K, and 
under an external field of 10 Tesla. 



APPENDIX C: PAIR-FLIP NUMBERS 

The number of pair-flips can be calculated with the 
pseudo-spin model, as an a posterior check of the pair- 
correlation approximation. The general formalism ia^ 

iVflip(t) ^ max^ |(T « E Kt ■ (CI) 

k k 

To be specific, let us consider the pair-fiip number un- 
der concatenated control. With the notations defined in 



Sec. VI the pair-flip number under the lih order concate- 
nation control is 



(0 

flip 



(-0 = E [(^ 



(C2) 



where {n^^ ^ and 



are the x- and ^/-component of 



the precession vector deflned in Eq. ( 69 1 , respectively, 
and T/ = 2't. With the condition t; ^ Bk ^ results in 
the leading order of the nuclear spin interaction strength 
are 









^ E^^ 




- E^ 

keGs 


in) 


- E^ 



r(0) 
'flip,. 



2 • 2 ^k 

.T smc — 
2 



(C3a) 
(C3b) 



(C3c) 

J2 Blrf sine" {Eut) < NBlrf. (C3d) 



k£l 



APPENDIX B: BOUNDARY OF MESOSCOPIC 
NUCLEAR SPIN BATH 



Intuitively, the nuclear spins in direct contact with the 
electron spin can be taken as the mesoscopic bath. The 
boundary of the bath is roughly defined by the condi- 
tion that the hyperfine interaction is stronger than the 
nuclear spin interaction for nuclei within the bath, and 
otherwise for those without the bath. The interaction 
between the nuclear spins within and without the bath 
across the boundary, however, could make this definition 
somehow arbitrary. One can define a larger spin bath by 
including one or more layers of nuclei outside the bound- 
ary. Since the nuclear spin interaction is much weaker 
than the hyperfine interaction, the cross-boundary inter- 
action should be unimportant in timescales of interest. 
Here we present the numerical test of the dependence of 
the electron spin decoherence on the choice of the bound- 
ary of the mesoscopic nuclear spin bath, particularly by 
examining the real-time behavior under the Carr-Purcell 
control. The result in Fig. [20] confirms the assumption 
that the slight ambiguity in defining the mesoscopic bath 
as a closed system is unimportant. 




t(|Lls) 



FIG. 20: (a) The electron spin coherence under the control of 
a Carr-Purcell pulse sequence, for various choices of boundary 
of the nuclear spin bath: For the dashed red, solid black, and 
dotted blue lines, the boundary is set 0, 1, and 3 layers of 
unit cells outside the boundary defined by the QD potential 
walls. The three lines are indistinguishable in the figure, (b) 
The deviation (amplified by 100 times) of the coherence for 
various choices of the bath boundary, referenced to the first 
choice of boundary. The QD is as in Fig. |19[ 
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FIG. 21: (a) The electron spin coherence in the Carr-Purcell 
control configuration, (b) The number of pair-flips due to 
the intrinsic nuclear spin interaction (dotted red line) and 
that due to the hyperfine-mediated interaction (solid blue fine, 
ampUfied by 10 times). The InAs QD is as in Fig. 19 with 
B = 2 Tesla. 
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FIG. 22: Schematics for the range of the QD size and the 
field strength (shadowed area) where the pair-correlation ap- 
proximation is justified in the timescales of interest in the 
single-pulse Hahn echo configuration. The upper and lower 
bounds set limits on the number of the intrinsic local pair-flips 
and on that of the hyperflne-mediated non-local pair-flips, re- 
spectively. 



Here the number of local pair-flips and that for hyperfine- 
mediated non-local pair-flips in FID conflguration are 
identified with the subscripts B and A, respectively. As 
illustrated in Fig. 11 (b), the pseudo-spin precession for 
non-local pair-fiips is reversed when the electron spin is 
fiipped, so under the concatenated control the number 
of hyperfine-mediated pair-fiips will oscillate periodically, 



with the maximum value ^ (t) at (2n-|-l)T and min- 
imum value at 2nr. At the spin echo time rj for Z > 0, 
the number of pair-fiips is contributed by the intrinsic 
nuclear spin interaction. The number of pair-fiips is plot- 
ted in Fig.[2T]for a typical QD in the Carr-Purcell control 
configuration, showing the features discussed above. 



To justify the pair-correlation approximation, the 
small pair-fiip number condition 



(C4) 

is to be fulfilled (see Appendix D]). Thus, the criteria for 
the timescales r and r; being short are given by 



N- 



NBtrf « 1, 

^n-^Alr^ « 1, 



(C5a) 
(C5b) 



for consideration of the intrinsic interaction and the 
hyperfine-mediated interaction, respectively. For an InAs 
QD containing about 10^ nuclei under an external field 
of i?oxt ^ 1 Tesla, the valid timescales is estimated to be 



Tl < 100 IIS, 

r < 10 ^s. 



(C6a) 
(C6b) 



A natural timescale of interest in experiments is the 
decay time of th e Ha hn echo signal, which can be esti- 
mated from Eq. ( 43b ) and ( 50 I to be 



H.B 



^5/125-l/2_^-l/2^ 



(C7) 



This timescale together with Eq. ( C5 1 sets ranges for the 



field strength and the QD size where the pair-correlation 
approximation is valid: 



N-'°/'n:'B^^Al « 1, 

as schematically shown in Fig. [22] 



(C8a) 
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APPENDIX D: ERROR ESTIMATION FOR 
PAIR-CORRELATION APPROXIMATION 

To estimate the error of the pair-correlation approxi- 
mation, we express the exact and the approximate nu- 
clear spin wavefunctions as 



exact, uncoiT/ 



\Jr 



± \ 
PCA,uncoir/ 



exact, corr/ ' 
PCA ,corr/ ' 



(Dla) 
(Dlb) 



respectively, where |Jetact,corr> and |J^cA,corr) denote 
the wavefunctions containing correlated pair-fiips, and 
l>7c^act,uncorr) and |J'pcA,uncorr) are the parts containing 
uncorrelated pair-fiips. The evolution starts with a ran- 
domly chosen nuclear spin state \J'). The uncorrelated 
part of the wavefunction in the exact solution and that in 
the approximation are determined by the same Green's 
function Go{t,t') and by the probability amplitudes at 
the initial state C^f^f^ j{t') and Cf^^^^ j{t') as 



\q± \ _ / ii\r<± 

-uncorr 



) = / Go(i,i')C^cA,j(t')|J)rfi'- (D2b) 
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After a few pair-flips have occurred while still well before timated to be 
the correlated part of the wavefunction becomes signif- 
icant, the probability amplitude of the initial state \J) 
decays already to zero. So the amplitudes at the ini- 5C\ _ ^ 

tial state are essentially the same in the exact and the 
approximate solutions. Consequently, the wavefunctions 
containing only uncorrelated pair-flips are equal in the 
two solutions 



I *-^xact , uncorr ) 



l>-^PCA,uncorr)- 



(D3a) 



(7- I 7+ 

X*-' exact, corr h-'exa 



J. 



PCA,corr 



7+ 

•-^PCAjCorr 



<Pr, 



(D4) 



Thus the error of the pair-correlation approximation in 
calculating the electron spin coherence correlation is es- 



where Pcorr is the probability of having correlated pair- 
flips. When the number of pair-flips iVflip <C N , Pcorr ~ 

1 — exp (^—qNI^^/N^ {q is a factor close to the number 

of nearest neighbors of a nuclear spin). 
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